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Summary 


A theory is proposed of the self "sustaining oscillations 
of a weak shock on an airfoil in steady, transonic flow. The 
interaction of the shock with the boundary layer on tho airfoil produces 
displacement thickness fluctuations which cojuvect downstream and generate 
sound by interaction with the trailing edge. A feed-back loop is 
established when this sound impinges on the shock wave, resulting in the 
production of further fluctuations in the displacement thickness. The 
details are worked out for an idealized mean boundary layer velocity 
profile, but strong support for the basic hypotheses of the theory is 
provided by a comparison with recent experiments involving the generation 
of acoustic 'tone bursts' by a supercritical airfoil section. 
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*!• Introduction OF POOR QUALITY 

The production of found by helicopter rotor bladee having 
transonic tip speeds is strongly influenced by the presence of shock 
waves near the blade tips (see, e.g. , refs. [l~3]). Acoustic waves 
are generated by (i) the interaction of the shock with inflow turbulence 
and/or the trailing vortices of the other rotor blades, and (ii) the 
intensification of the nonlinear quadrupole volume sources which occur 
in the Lighthill theory of aerodynamic sound [4]. In this paper we 
discuss a possible additional mechanism associated with an instability 
of the shock caused by its interaction with the boundary layer on the 
blade. This interaction is strongly dependent on the Reynolds number 
when the flow is transonic [5—7^ , the position of the shock being 
extremely sensitive to temporal variations in the properties of the 
boundary layer. 

Sued et al have reported [s] the occurence of intense, high 
frequency and highly directional tone bursts during acoustic tests of a 
scale model of a general aviation propeller operating at high subsonic tip 
speeds. The amplitude of the burst was of the same order as the propeller 
noise, and observation strongly suggested that the tones were produced by 
oscillating shock waves on the blades. The phenomenon was absent if flow 
separation occured at the blade tips and if the boundary layers ahead of 
the shocks were made turbulent; this is presumably a particular manifesta- 
tion of transonic buffeting [9] . The extreme sensitivity of the shock to 
perturbations in the ambient flow is also illustrated by the experiments 
of Tijdeman [10] using an airfoil with an oscillating flap. These 
demonstrate that under certain conditions the shock leaves the airfoil 
and propogates as an intense acoustic wave to the far field, an effect 
which is in qualitative agreement with theoretical work of Williams JjLlJ . 
Similarly, Magnus and Yoshihara Q.2] argue that discrepancies between 


1.1 


experiment and their numerical prediction* of motion cau*ed by a pitching 
airfoil in transonic flow are probably a revolt of the atrong interaction of 
a ahock wave with the boundary layer. 

An idealized transonic flow problem i* investigated below in 
order to examine a mechanism of ahock wave/boundary layer interaction 
which i* possibly responsible for the intense acoustic fields observed by 
Succi et al [8], The analysis proceeds from the unsteady, transonic flow 
equations which are linearized about a steady mean flow with account taken 
of the displacement of the shock, in the manner described by Williams [ll]. 
The model problem is illustrated schematically in Figure 1. A two- 
dimensional flat-plate airfoil is placed at zero angle of attack to a mean 
flow. Above the airfoil the mean flow is uniform and supersonic ahead of 
a weak shock wave which in the undisturbed state ri assumed to be normal 
to the airfoil and of unlimited extent. Time harmonic oscillations of 
the shock about its mean position are considered At the root of the 
shock the motion induces fluctuations in the ^ispia cement thickness of the 
boundary layer which propagate downstream as a surface wave on the airfoil. 
At the trailing edge this wave feeds into an unsteady wake and it 

responsible for the production of edge-generated sound which subsequently 
interacts with the shock thereby inducing further fluctuations in the 
boundary layer. Self-sustaining oscillations are possible provided the 
returning sound waves are in an appropriate ph&se relation with the motion 
of the shock, and are of sufficient amplitude. These conditions determine 
a preferred wavenumber for the bovmdary layer waves, and our objective is 
to determine the frequency of the oscillations in terms of the characteristi 
of the mean boundary layer. 

The analytical problem has similarities with that investigated 
recently by Goldstein et al Q.3] , concerning the instability of ahocka 
of arbitrary strength in cascades, although no account was taken of 
shock wave /boundary layer 


is assumed to be sufficiently weak that the generation of vorticity and 
entropy by its motion can be neglected. The unsteady flow in the main 
stream behind the shock may accordingly be expressed in terms of a 
velocity potential 4 . At the outer edge of the boundary layer the 
normal derivative of 4 equals the displacement velocity of the 
unsteady boundary layer. Application of this condition is simplified 
by wans of an approximation analogous to that used in thin airfoil 
theory [l4] , which permits it to be imposed on the upper surface of the 
airfoil (c.f. [15-17]), 

The principal difficulty in formulating the problem is the 
modelling of the production of the displacement thickness waves. The 
approach adopted here is based on the hypothesis that the functional 
form for the transition in the boudary layer structure across the shock 
is invariant in a frame of reference fixed with respect to the root of the 
shock. This assumption of quasi-static behaviour at the root permits the 
unsteady motion in the boundary layer to be determined in a linearized 
approximation provided the mean properties of the boundary layer are known. 
Results are given here only for the highly simplified case in which the 
mean 'velocity profile is approximated by a step function, a procedure which 
has been successfully exploited by Ffowca Williams and Purshouse [l 8 ] and 
Goldstein [l9] in analytical studies involving unsteady boundary layers. 
Evidently this approach is also applicable to other shock vave/boundary 
layer interaction problems. For example, it provides an excellent starting 
poire for investigating the interaction of blade-tip generated shocks of a 
ducted rotor with the boundary layer on the walls of the duct. 

The basis of Che model problem is discussed, formulated analytically and 
solved in i§2,3. In §4 the equation for the characteristic frequency of 
the self-sustaining oscillations is obtained; numerical results given in 
f5 are examined in relation to the experiment of Succi et al [b] . 
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12. rpnnulation of th« aodel problem OF POOR QUALITY 

A two-dimensional rigid airfoil occupies the portion 
- c < Xj< 0 of the Xj-axia of a rectangular coordinate uyatem (x^Xg), 
in the presence of a Mean flow in the positive x^'diraction. In the 
undisturbed state a normal shock wave extends from the upper surface of 
the airfoil at x^ ■ HI to x^ - (Figure 2(a)). Attention is confined 
to the veak-shock/transonic regime in which the speed U of the main flow 
downstream of the shock is constant and 

fsM-l-|! 2 / 2 , (l) 

where c is the speed of sound, and B ^ /l-M 2 ^ « 1 . The Rankine- 
Hugoniot relations applied to a weak shock imply that the main flow 
Mach number M_, say, upstream of the shock is given by 

H.-1 + B 2 / 2 . (2) 

Motion of the shock wave produces fluctuations in the boundary 
layer and wake which in turn react back on the shock. It is assumed that 
all such back-reactions occur via acoustic path6 intersecting the shock 
from behind, i.e., possible interaction channels involving with the passage 
of Sound around the leading edge of the airfoil, after propagation through 
the subsonic flow below the airfoil, are ignored. Similarly, the trans- 

mission of acoustic disturbance into the supersonic region via the subsonic 
boundary layer is neglected. 

Consider a time-harmonic perturbation proportional to e . 

Let the instantaneous position of the shock be represented by 

Xj - - l * x(x 2 )e” 1, * )t , (3) 

where *(x 2 ) is assumed to be sufficiently small that the perturbation 
equations may be linearized. For a weak shock the motion in the main 
flow above the airfoil is described by a velocity potential $ e , 
in terms of which the perturbation velocity v « Vf (the exponential time 
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factor it here and hcncaforth auppraaaad) 


In a linearised approximation 


4 satisfies 
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- ik o*“£ 


a 2 


a 2 


a*i 2 ?x 2 2 


<4) 


* In the subsonic region > - i ♦ , where k o - u/c . 

In the supersonic flow ahead of the shock the motion is steady, and 




| 

I 




4-0 (x 1 < - l 4 e(x 2 )) (5) 

The potential must also satisfy the following conditions at 
the undisturbed location x^ - - l of the shock: 


li ♦ 0 . 

ax ’ 


ue* 


*(x«) - — ^ , 


-9 

cp- 


( 6 ) 

(7) 


where it is understood that these conditions are to be satisfied by 4 and 
S4/3x^ as x 1 ■+ - i from the downstream region. The derivations of 
these formulae are given by Williams [ll,20]. 

In the (shock free) main flow region below the airfoil the 
motion is assumed to be subsonic everywhere and 4 is required to satisfy 
equation (4) . 


Boundary conditions at the airfoil and wake 
Let v q (x^) denote the normal (i.e. x 2 -) component of velocity 
on the surface x 2 - 6 lying just outside the boundary layer on the upper 
surface of the airfoil, 6 being the boundary layer thickness; this will 
be referred to as the boundary layer displacement velocity . On x 2 - 6 
4 must satisfy 

‘ V n <x l > ' (I 1 * °> • (B) 


2.2 
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In order to snake the subsequent analysis tractable, this condition 
ie imposed on x 2 - * 0 , an approximation which is impacted to be 
valid provided the length acalc U/w of boundary layer disturbances 
is large relative to 6. Boundary layer displacement fluctuations are 
assumed to be absent on the lower surface (x 2 • - 0) of the airfoil, 
where it is accordingly required that 

|i « 0 , - C < % x < 0 , * 2 • - 0 . (9) 

He shall actually apply this condition over the semi-infinite interval 
- • < < 0 . this will ensure that acoustic disturbances cannot 

impinge on the shock from the supersonic region, and avoids difficulties 
(which cannot be incorporated i&to the present idealized model) arising 
from the fact that the mean flow ahead of the shock must actually vary 
with position* 

Similarly, if tdh/U « 1 , where h characterises the mean 
thickness of the wake downstream of the trailing edge, the perturbation 
pressure 

p - p 0 (iw - U — •) * , (10) 

(p 0 being the mean density, taken to be constant), may be assumed to be 
continuous across the wake. In the usual approximation of thin airfoil 
theory [l4j , this condition may be imposed on the centre-line • 0, 

Xj > 0 of the wake, and implies that 

r i iwx l /U 

U] “ Ac (% x > 0) , (11) 

where £f] ■ $(x^, ♦ 0) - ♦(x^, - 0) defines the, discontinuity in f 
across the x^-axis. The value of the constant A is a measure of the 
strength of vorticity shed from the trailing edge, and is determined by 
the requirement that the perturbation pressure and velocity should remain 
finite at the edge (Kutta condition) . 






Q* 


A non-xero value of • A l» associated with the presence of 
an asymmetric (sinuous) disturbance in the wake. Symmetric "breathing" 
modes are als*» possible, however. When wh/U is small the pressure 
is continuous for such modes,, but the normal velocity exhibits e simple 
discontinuity across x^ * 0, > 0 . In this case the phase velocity 

1. aqual to th. .doit,™ vlocity «jr, in th. **k. [2l] , mi 

if tc s » U^/w , the breathing mode satisfies 

iK .*i 


’%] ‘ «* 


<X l > 0) • 


( 12 ) 


where B is a constant. The value of B is obtained by requiring that 
there be no net flux of fluid from the wake and boundary layer. This 
condition is derived from the equation of continuity on the basis that the 
characteristic acoustic wavelength is large relative to the thickness of 
both the wake and the boundary layer. Thus in the absence of displacement 
velocity fluctuations on the lower surface of the airfoil, B is related 
to the displacement velocity v r by 


f 

-C 


v n<‘l )dx l 


j* Be 


1K . X 1 


dx, ■ 0 


( 13 ) 


Convergence of the second of these integrals is assured by assigning to 
u) * Uk # a small positive imaginary component which is subsequently allowed 
to vanish. The validity of this procedure is a consequence of the causality 
condition which requires all field variables to be regular for sufficiently 
large and positive Im u . The disturbances in the wake arise as a result 
of the displacement velocity fluctuations in the boundary layer and of the 
motion of the shock. Hence 

B - ix # j° v n (x^)dx^ . (14? 

-C 
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Modelling the interaction of the .hock end boundary ” QUALITY 

layer 

Oscillatory notion of the shock lit frequency u generates 

displace nent velocity fluctuations whose specification according to linear 

theory depends on a knowledge of the nean boundary layer velocity profile. 

Ve consider here only the highly idealised velocity profile illustrated in 

Figure 2(b), in which the nean shear in the outer region of the 

boundary layer is concentrated into a vortex sheet at a stand-off 

distance 6 from the airfoil. If the flow in the inner region 

0 < x 2 < 6 has speed V in the x^-direction, the vortex sheet can 

A i« ± x. 

support displacement velocity waves proportional to a , where 



pfc'&vided w6/V is small, i.e., that the wavelength is large relative to 6 
(c.f. ref. [l6] , 53) . 

In the undisturbed state there will exist a length scale d 
determining the distance in the x^-direction over which the mean properties 
of the boundary layer change across the shock. When oscillations occur 
at frequency u it is anticipated that at distances exceeding d downstream 
of the shock the unsteady motion in the boundary layer consists of a linear 
combination of the two displacement velocity waves determined by the 
dispersion equation (15). This notion leads to the quasi-static 
representation of the interaction of the boundary layer and shock described 
below. 

In the steady state the boundary layer thickness fi increases 
rapidly with x across the shock, and we write 


6 - ^(xj* l) 


(16) 
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and uium that 

4 

rCitj+ft) %% for A < - d 

% i ♦ «q /5 o for V 4 * 4 

* t0 that 6 - 6 0 , « 0 ♦ *0 respectively in the vegions upstream and 
downstream of tha shock. Tha generally much slower variation* of fig*®]) 
with x dua to tha natural growth of tha boundary Idjrer ara ignored. 


In generalising equation (16) to unsteady action, observe that In 
practice ’’fanning" of the shock aust occur in the vicinity of the 
boundary layer, and therefore that the shock- location aquation (3) 
becomes ill-defined as Xj * 0 • This fanning will be neglected 
to the extent that when x^ 0 it will be assumed that a representative 
value of Sq ■ »(0) can be defined to determine the instantaneous 
position of the root of the shock just outside the boundary layer. 

Let c(x^,t) denote the unsteady boundary layer thickness, 
and introduce the representation 


C « 


F(x+l) - * Q 


r 


ik(x,*H) 

(k-XjHk-xj) 


dkV , 


(18) 


where k^, k* respectively denote the wavenumber k defined by equation 
(13) when, the -, ♦ sign is taken. The path of integration in (18) runs 
below both of the poles at k - Kj, k* (this can always be ensured by 
taking w to lie in a suitable region of the upper half-plane). The 
function f(kd) is assumed to be independent of u , to be regular on 
the real k-axis, and to vary significantly only for variations in kd 
which are 0(1), . In practice d and 6 Q are of cotq>arable magnitudes. 


r 
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mad both arc hero required to be small relative to the wavelength 
<vl/Re(ic, ? ) of the displacement velocity waves, This iaplies that the 
singularities of f(kd) in the complex plane lie further from the real 
axis than the simple poles at k * K x ,K l * Accordingly, when 
x^ 4 l >> d > 0 the value of the integral in (18) is dominated by the 
residue contributions from these poles, and the notion thereby described 
is entirely associated with the corresponding displacement velocity waves. 

When 4i> is real, respectively determine boundary layer disturbances 

which grow and decay exponentially with increasing x^ • 

In the opposite extreme in which |Xj+i] < 0(d) , the 
displacement velocity poles make an insignificant contribution to the 
integral, and (18) reduces to 


i 


C 


mx+o - * 0 


r 


ik (x.+l) 

f (kd)e 1 dk - . 

, 


( 19 ) 


The hypothesis that the structure of the boundary layer in the neighbourhood 
of the ahock is invariant in a reference frame fixed relative to the shock 
requires that f(kd) be interpreted as the Fourier transform of 


F'(Xj) ■ dF/dx, i.e.. 


■ h ^ • 


f(kd) 


in which case (19) becomes 


C • 6, 


F(x 1 <’l-* () ) ♦ 0 (eg) 


( 20 ) 


( 21 ) 


In the following discussion the functional form of F(x+1) is assumed to 
be known. 


2 . 7 
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The displacement velocity v (O' may now be expressed In terns 

u I 

of the displacement, * 0 of the root of the shock by substituting the 
representation (16) into the linear theory foraula 



C 


<*! > - D 


( 22 ) 


giving 


W ■ 1 Vo 


r 


ik(x.+l) 

k 2 (urUk)f(kd)e 1 
(k-KjXk-icj) 


dk 


(23) 


Similarly, since F f (x) vanishes for x^ ♦ 1 < - d , upstream of 
the shock, equation (14) becomes 


1 Vo*. 


ik£ 

P kf(kd)e 

J Tk^)Tk-K*y 


<Sk . 


(24) 


Boundary value problem for d 

Collecting together the principal strands of the above discussion, 
the boundary value problem for the potential + may be stated thus: 


For (Xj > - 1, x 2 > 0) and (- •• < x 1 < •, x^ < 0) , 
admissible solutions # of the equation 


find 


(-“.•» si:)' -fe- 


iX 2 


( 25 ) 


which satisfy the radiation condition at large distances frsa the airfoil 
together with the following: 


I* Conditions on the Airfoil: 
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■ i * 0 4 0 f 


ik(x.+£) 

- , i k 2 (w-Uk)f(kd)e 

Vo j * • 


(- t < Xj « 0, *2 - ♦ 0); 



(- • < Xj < 0, Xj • - 0) . 


II, Condition! in the vake (x^ ■ 0, x^ > 0) ; 


W 


iwx. /U 

A e 



. 

, 


III. Conditions at the shock: 

*0 “ <Jt i *J ’ 0) • 

The coefficient A is to be deterained by application of the Kutta 
condition at the trailing adge of the s>foil, and B is given in terns 
of Zq by equation (24). 

When used in conjunction with the dispersion equation (15), the 
aolution of this boundary value problen provides the characteristic 
aquation for the admissible frequencies u for which self-sustaining 
oscillations of the shock wave are possible. 


(26a) 

(26b) 

(27a, b) 

(28a) 

(28b) 


S3. Solution of the boundary value problem 
The Fraadtl-Clauert transformation 

l KMX 


* 

X 


♦ •' 

k O /0 


X ■ Xj/8 
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( 29 ) 


reduces equation (25) to the standard form 


ii ♦ 

3X 2 3x 2 

2 


♦ - 0 , 


(30) 


sad the boundary conditions (26) - (26) become: 


X. Conditions on the airfoil: 


W * j- * r v 2 (urUv)f(vd)e i ^ pg,> ' v t ^ 
3x^ 1 0°0 j (v-Kj) (v-icj) 


ou 


(- L < X < 0, x 2 - ♦ 0) . 


where L - 1/8 , o - v8 ♦ KM • 


(31a) 


34) 

3x, 


(- • < X < 0, * 2 • - 0) . 


(31b) 


II. Conditions in the wake (x,j ■ 0, X > 0) : 


M - A e iKX , (*c - u/BV) i 


W iic X _ 

-Be , (k - k 6 ♦ KM) 

■ s 


(32a) 

(32b) 
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III. Conditions at th« shock: 
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|i ♦ iic* “ 0 (X -► - L, *, > 0) { (33a) 

* - - — * (X * - L, 0). (33b) 

0 cB* 2 

Equation (33a) has been ainplified by asking use of the condition 

0 2 « 1 . 

Since the problem is linear we can set 


♦ • *i ♦ * 2 » 

where ^ represents the flow induced by the boundary layer displacement 
velocity fluctuations, and i> 2 i* that produced by the oscillating 
motion of the shock. The functions * L , i> 2 BU,t reduced forms 

of conditions (31), (32), and their sum, ij/, must satisfy (33). 


Bounds x/ value problem for ^ 

The perturbation due to the displacement velocity fluctuations 
is determined by ignoring the presence of the shock, i.e., by assuming 
the flow to be subsonic everywhere and by imposing (31a) over the semi- 
infinite interval (- « < X < 0) . This procedure is permissible 
provided that the conposite potential satisfies conditions (33) on 
the shock. 


A Fourier integral representation for is obtained in Appendix 


♦l" 


-Hopf technique [ 22 ] . 

In 

* 2 > 0 

Vo f'l 

r v 2 (w-Uv)f (vd) 

, % 

► 

► 

1 

47T j J 

— • 

(v-O (v-<J) 


“TTiOi 

» • 

K /V 

6 

L i 1 

* 

l 


k - a + lO 


fk-^-iO) 




k-o+iO k-K-iO 


i{vl+kXn(k)x 2 > 


x e 


dkdv , 


(35) 


in which the notation t iO indicate* that the contour of integration 
in the k-plane passes above the pole at k * a and below the pole* at 

k * Sc • k • The function yOO i* defined by 
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where branch cut* for the radical* on the right-hand aide are taken 
respectively in the upper and lower halve* of the k-plane w such that 
when K ia real and positive, y(k) is positive on the real axis for 
}k| < K and positive imaginary for |k| > K . 

Boundary value problem for 

Motion of the shock causes sound to be radiated in the downstream 
direction. Since the influence of the displacement velocity is included 
in • Che shock associated sound must satisfy 


3*2 

3X 2 


0 


(- - < X < 0 


0 ) 


(37) 


rather than (31) , where again the condition is imposed on the half-line 
(- “ < X < 0) . In the wake there is no need to account for breathing 
modes generated by the displacement effect upstream of the edge, and ^ 
is therefore required to satisfy conditions (32) with B - 0 . 

In the absence of the edge an appropriate representation of the 
motion in Xj > 0 would be 


|} 

V 

R 


*2 | T(X)cos(Xx 2 ) e iya)(X+L) dX 

0 


» 

(38) 


If 


*2 



(39) 


is the diffracted field produced by the edge and wake, and may again 
be determined by the Wiener-Hopf procedure. In Appendix 2 it is shown 
that, for x^ > 0 

s^~ 

4 


3.3 



miGiNM- Pf (*■ 
■jp poOR Q uAUT 



k=7TxRo 


| i(kX+Y(k)x 2 4 ir(X) L >^ 


l 


( 40 ) 


Equations (38) -(40) data mine ^ in terns of the as yet unknown 
function ?(*) . 


The integral equation for Y(A) 

Condition (33a), to be imposed at the undisturbed location of 
the shock, provides an integral equation which determines Y(A) in terms 
of the displacement z^ of the root of the shock. To obtain this 
equation we introduce the Fourier cosine transform I(y) of a (function 
gCx^) defined by the reciprocal relations 


I(p) * j 
0 

g(x 2 ) - \ 


g(x 2 )cos (yx 2 )dx 2 

J" g(p)cos(yx 2 )dy 
0 


t 


(41a, b) 


(see, e.g. Erdelyi et al J23j) . 

Substituting ^ ♦ t|> 2 into (33a), using the representations 

(35), (36), (40), and taking the cosine transform, we find that Y 
satisfies 


f(M) * 


i 

2* 


x 


f*-Y(n)| r 

UnfoJ J 

o 

(y(A)*hf(wi 


W Y(y) “ 


1 


e iL(y(u)n(^)) dX 



t 


1 

(k+y ( u) ) 


X - - L 


(42) 


where Che cosine transform is given by 



1 Vo 


r v 2 (<r-Uv)f (vd)e^ Vt *°*^dv 

j — 


+ 


id a 
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v 2 (srUv)f (yd) 
dv-Kj) ( v-kJ)y ( y) 
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Y(|i)-H7 


« 

. tjnTu? i i 

/jjT ^ — \ |y(m) + k “ ~v(\i)+c 


l e i(vi-Y(w)X) dv 


(43) 


The various terms in equations (42), (43) have the following 
interpretation: 

Y(y) is the contribution to (33a) from Yj » th* direct 
field generated by the shock; the integral involving Y(X) on the 
left of (42) represents the contribution Y^ of sound waves generated 
at the trailing edge by diffraction of the primary fiel d j’T, • The 
first term on the right of equation (43) denotes the influence of the 
boundary layer displacement velocity fluctuations when the edge and 
wake are ignored; the second term represents the effect of sound 
generated at the edge due to its interaction with the displacement 
velocity wave* 


14 » The characteristic frequencies of self-sustained oscillations 


A second relation between * Q and f(X) can be obtained from 
condition (33b) imposed at the root of the shock, vis: 
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t iKML 

cP 2 


^(-L.O) *► l£(-L,0) ♦ T®(HL,0)} 


(A4) 


The elimination of * Q , T(X) between (42), (44) provides tax equation 
relating the frequency m and the conjugate displacement velocity 
wavenumber* • Tbs roots of this equation define those admissible 

values of Kj,** for which self-sustained oscillations are possible, and 
when substituted into the dispersion equation (15) yield the characteristic 
equation determining w in terms of the boundary layer thickness 
6 ■ 6 q + 6q downstream of the shock. 

Self-sustained oscillations are envisaged to occur as a result 
of the following feed-back loop of mechanical processes (see Figure 3) : 

1. The displacement Zq of the root of the shock generates 
displacement velocity waves in the boundary layer; 

2. Sound waves are produced by the subsequent interaction of 
the displacement velocity waves with the trailing edge; 

3. The impingement of the sound on the shock closes the 
loop by generating further disturbances in the boundary 
layer . 

These observations lead to the following simplifications of equations 
(42), (43) and (44): 

In equation (42) the second term on the left hand side is 
discarded. This expresses the influence on the shock of sound waves 
orignally generated by the shock (^) and subsequently diffracted at 
the trailing edge. The amplitude of these waves is expected to be small 


4.1 


'i * 


coexisted with the sound produced by the boundary layer/edge interaction 
because of the exponential growth of the displacement velocity waves 
as they propagate towards the edge. 


. vlt 






•V 


When the second integral on the right of equation (43) is 
i inserted into the right hand side of (42), it describes the effect on 
the shock of sound generated by the boundary layer/edge interaction. 
This interaction will be dominated by the exponentially growing 
displacement velocity wave (proportional to e I 1; , whose effect is 
determined by the residue contribution to the integral from the pole 
at v " . This will be taken to represent the principal component 

of the edge generated sound. 

Taking account of these approximations we accordingly reduce 
equation (42) is reduced to 

*(v0 - z 0 <x 1 (v) ♦ x 2 (vO) » 


where x 1 X 2 » respectively denote the local and edge-diffracted 
influences on the shock of the boundary layer motion, and have the 
explicit representations: 


• ~ J.CV4XLI . 

1 Q C r v 2 ((p-Uv) (o+K)f (vd)dv 

<V“< i )(v-x*)(K 2 -u 2 -o 2 ) 


ft 

i 



where 


x 2 (v) 


Trg O fx“Y(u)' 

2 \K+r(v y ( 


2 i(< l+y(y)L) 

Xj(ur-UK I )f(K I d)e 


(*!-«*) 


• 

. k s /k i i t I'E+Ttu)' 

f 

1 1 

y(?j)+K s Y(V) +0 I 

|y(v)+< Y(y)-o I j 

. 




(45) 


(46) 


i (47) 


( 48 ) 


4.2 


Turn attention now to equation (44) * The first term in the 


brace brackets on the right hand aide accounts for the backreaction 
of the boundary layer on the motion at the root of the »hock. Aa above, 
this decomposes into two components, the first describing the local influence 
of »ihe boundary layer (when edge effects are ignored), and the second the 
effect of acoustic waves produced by the boundary layer /edge interaction. 

Both effects are contained in the double integral representation (35). 

When X *♦ - L on " + 0 the integration contour in the k-plane may 
be displaced to - i« . In so doing the integral along the real axis is 
transformed into the sum of two terms Z Q®2 * ,ay * respectively 

equal to the residue of the simple pole at k - o (which characterises 
the local effect of the continuum of boundary layer displacement velocity waves 
specified by equation (23)) and an integral around $ contour enclosing the 
branch cut of , which extends from k " - K to - K - i" , The 

latter represents the edge generated sound. 

Thus we find: 


« _ r _-iKML f° v 2 (w-Uv)f (vd)dv 

®1 V J T ( 0 )(v-« ) • 

•HB 


(49) 


at J the branch-cut integral can be reduced to the form 


■ f d T 


J(v f Oe“ KU 

7 F~ 


d C 


(50) 


where ^(v^O is regular on the £-ax is and 'v 0(1/0 «s 5 + ® . 

Now B 2 << 1 , so that KL 5 k^i/B 2 is a large parameter which 
characterises the foreshortening of the wavelengths of sound waves generated 
at the trailing edge and propagating upstream against a near sonic mean flow. 
To leading order, we therefore have 


0 2 * f *< v »°> dv f jz d * 

mm 0 ' 
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Jkf 


0) dv 


( 51 ) 


Tb« remaining integral with raspect to v is approximated by the residua 
contribution from the pole at V * Kj (see equation (35)) on the basis , 
discussed above, that the dominant edge-generated sound is produced by 
the exponentially growing displacement velocity wave proportional to 
e l 1 • In this way equation (51) is found to have the explicit form: 


( _ k2(w-Uk t ) 

e 2 * ~ MEuJ 1^-kJ) f(K i d) ■ 

1 , ^ ( 1 1 1 
I* 4 * " 


"s /K I 

(K+K g ) 


i [kjI+KL-it/4] 


(52) 


This and equation (49) determine 


V‘ l -°> ■ "o{v*k} • 


(53) 


The second term in the brace brackets of equation (44) is the 
contribution of the primary field due to the notice of the shock, and 
it follows immediately from (38) that 

- f f *(X)dX (54) 

0 

The final term in (44) is discarded for reasons discussed 
previously, since it represents the influence of sound produced by the 
diffraction of ^ at thc trailing edge. Hence substituting from (53), 
(54), equation (44) may finally be cast into the approximate form 
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When this expression for t Q la inaertad Into aquation (45), 
the following integral equation ia obtained: 


2. lrML /»cBM(x 1 (w)+X 2 (>'>) r 
.rm, — J 

[l ♦ (G^) 0 


f(u) ♦ 


¥(X)dX - 0 , 


(0 < y < •) . 


The characteristic equation relating admissible values of u, , k* 
ia now deduced by integrating this result over 0 < y < •• : 


iKML 


' I (V*z) + f ( *1 
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(y) ♦ x 2 ^)) d w 


| « o . 


Define 


A more convenient form of this equation ia obtained as follows. 




iKML f n$ 


,2 2 


f x^yjdyj 
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v 2 (m-Uv)f (vd)dv 


[c) (v-KjHv-kJ, 


- 1 r dM r — 

0 Lm (v-Kj) (v-K* 


(ur-Uv) (o+Qf (vd)dv 
( v - k *) (K 2 -y 2 -o 2 ) (icn(y)) 


where use has been made of equations (46) , (49) . 

The presence of the factor e^(y)L £ n rftpresen tation 
(47) of x 9 Cy) indicates that, when KL - i/6 2 is large, as assumed 


e, the value 


t (y)dvi may be approximated by the method of 
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stationary phase [ 22 J • Taking account of this remark ve define 




i8 ^ 2e 


iKML 


vc6‘ 


* G 2 . 


X 2 (y)dy 


where (from (47) anr* (52)) in a leading approxisuition 


and 



< 59 ) 


(60a) 


(60b) 


The characteristic equation (57) say accordingly be set in the 


form 


1 * A (iCjjKjpO)) ♦ B (tj.icJ.w) 


ie 


(61) 
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15. Discussion of the Characteristic Equation OF POOR QUALITY 

The characteristic aquation (61) gives the condition which 
wist be fulfilled if self-sustaining oscillations are to occur at 
frequency u for given values of the displacenent velocity wavenumbers 
K I* K { * 7)16 f * r,t two tena » oo the lefthand side represent the 

influence of the motion at the root of the shock when diffraction at 
the trailing edge is ignored. The term in 8 gives the backreaction 
at the root of sound produced during the passage of the unstable 
boundary layer wave over the trailing edge. This interpretation is 
evident from the structure of 

k l 

Be 0 - ReCKjl) + ~ “ , (62) 

wherein (i) Rc(kjI) is the phase change associated with the convection 
of a boundary layer disturbance from the root of the shock to the 
trailing edge, and (ii) k ft i/l-M) - tr/4 is the phase change experienced 
by sound waves radiating upstream from the edge to the root together with 
a correction (- ir/4) due to the cylindrical spreading of the waves. 

In order to derive quantitative predictions from equation (61) 
it. is necessary to introduce an explicit representation of the function 
f(vd) given by equation (20). Consider first the value of f(Kjd) 
which occurs in the definition (60a) of B . By hypothesis, both the 
transition length d and the boundary layer thickness 6 m 6^ * 6^ 
are small relative to the wavelengths of the displacement velocity waves. 
Since F'(x^) is significantly different from zero only for JxjJ < 0(d), 
it follows from (17) md (20) that, to leading order, 

s** 6 ' 

f(Kjd) s 2 iT J F , (» 1 )dx 1 * , (63) 

and this estimate will be e alloyed in calculating 8 . 


5.1 


To estimate the corresponding values of the coefficient A 


from its definition (58) , it will be assumed for simplicity that 
F'Cx^+t) is an even function with respect to the undisturbed location 
* - l of the root of the shock. When K^d « 1 the leading order 
s asymptotic approximation to equation (58) reduces to 


A (x 1 # »cj,u) 


r 



( 64 ) 


Numerical results are presented below for the case 


f(A) 



2w(H-A 2 ) 


» 


for which, 


A (x-.xi.u) 




This and ather expressions appearing in the characteristic 
equation (61) will be expressed in non-dimensional form in terms of 
the reduced frequency S and wavenumber W defined by 


(65) 


( 66 ) 


S « k t = — 
o c 

*1 

w = U ♦ iw - — 

1 lW I k • 


(67a, b) 


where are the real and imaginary parts of W . If the mean 

velocity in the wake is assumed to relax rapidly to that of the mean stream 
U , the wavenumber of the breathing mode in the wake is equal to 

w/U , and this value is used below. 

Equation (61) may now be expressed in the form 
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- j ♦ arg G 
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(68a, b) 


where 
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1 M 2 

(w-l) 

r mj 

r 

1 w iw/T 

M _ ] 

L J 

(i+m) 

(1+M) (1+M+B 2 W) [gZw-UM]* 

lfM (i+M+6 2 Kl| 



and n ia an integer. Equation (66a) ia the feedback loop condition that 
the total phaae change around the loop ahould be a multiple of 2ir . 

The term arg (G) ia the change in phaae introduced into the loop 
during the production of diaplacement velocity and aound vavea ; the 

remaining terms on the left of (66c) account for the effects of propagation 
discussed above. Equation (68b) determines the growth rate of the 
instability wave in the boundary layer , v* ich muat be auch as to ensure 
that the amplitudes of the sound vavea returning to the shock are sufficiently 
great to sustain the loop. 

Equations (68) are to be solved simultaneously with equation (15) 
relating the displacement velocity wavenumber to the characteristics of the 
assumed boundary layer velocity profile. In terms of dimensionless 
variables equation (13) gives for the exponentially growing wave 


t 





1 

v^EWu-v) 2 ] 1 ’ 


(70) 


where 

H m v/U 

is the ratio of the inm?r velocity of the model boundary layer (Figure 2) 
to the main flow velocity, and 6 ■ is the boundary layer thickness. 

Equation (70) characterises a disturbance which grows exponentially as it 


propogates at phase velocity V downstream fro* the shock., 
Numerical Results 
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Equations (68) , (70) are used to determine the dependence 

of the reduced frequency S of the self-sustained oscillations on the 

boundary layer thickness ratio 6/1 . The calculation is performed 

for fixed values of 6'/** d/6 * , H and V . The values of $ and 

o o 

6/1 determined by this procedure are actually discrete* and parameterised by 

the integer n (S increasing with n) . Their variations with n turns out to be 

sufficiently smooth* however, for it to be convenient to present the 

results graphically as continuous plots in the (S, 6 /l) -plane. The 

results illustrated in Figures 4,5 correspond to the case 6^/1 - 0.0005 

and M • 0.99. The range of variation of 6/1 is set by the constraints 

(i) that 6 * 6 o >o,the boundary layer thickness upstream of the shock* 

and (ii) the upper limits of validity of the approximate dispersion 

equation (70) (wherein the second term in the brace brackets is required 

to be small relative to unity). The solid and dashed curves are for 

d/6^ ■ 1 and 10 respectively; V is equal to 0.2 in Figure 4 and 0.6 

in Figure 5. The lower of these velocities is intended to model the 

order of magnitude of the phase velocity of the long wavelength 

disturbances in the boundary layer [24]. Prominent, large amplitude 

boundary layer structures tend to convect at about 60Z of the main flow 

velocity, and this is the situation modelled in Figure 5. In both cases 

it is evident that the predicted variation of S with 6/1 is not 

significantly dependent on the transition width d/6\ 

In the experiments reported by Burdges [7] the increase in the 
displacement thickness 6* of the boundary layer across the shock is 
typically of order 0.0011 . If, follcwiag Ffowcs Williams and Purshouse 
M . we identify the boundary layer thickness 6 of the vortex sheet 
model with the displacement thickness of the real boundary layer* 
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i.e., we take 6^/1 ■ 0.001, the results depicted in Figure 6 are obtained 
when V - 0.2. It is of interest to consider this case in relation to 
the experimental findings of Succi et al Ijfj . In their experiment 
A i 0.8 cms and acoustic tone bursts were observed at a blade referenced 
frequency of 12.6 KHz. This is consistent with the present theory in as 
much as in Figure 6 the resonance frequency ■ Sc/2*t varies between 
10.2 - 21.6 KHz when d/6^ ■ 1 , and between 8.1 - 18.3 KHz at the 
higher value d/6^ - 10 . 


56. Conclusion 


A linearized theoretical model of the interaction of a normal 
•hock wave with the boundary layer on a supercritical airfoil has been 
discussed. The system can oscillate at certain discrete resonance fre- 
quencies provided an appropriate feedback** loop condition is fulfilled. 

This is dependent on the convection velocity and growth rate of dist- 
urbances generated in the boundary layer by motion of the shock and on 
the subsequent interaction of those disturbances with the trailing edge 
of the airfoil. Resonance frequency predictions are consistent with 
measured "tone burst" frequencies observed in the study of the noise 
produced by a transonic propeller. 

The analysis of the coupling between fhe shock and the boundary 
layer is based on the hypothesis that the motion in the neighbourhood of 
the root of the shock is steady in a reference frame moving with the root. 
There, are obviously many differences between the model and real flows: 
shocks are not usually weak nor of infinite extent, are often accompanied 
by a separation bubble and never enter the boundary layer in one front. 

In addition, disturbances generated downstream of the shock leak through 
the boundary layer and can also propagate around the leading edge of the 
airfoil to modify conditions ahead of the shock. However, the theory 
provides a useful first 6tep in the understanding of unsteady shock/ 
boundary layer interactions, and would also be applicable to other aero- 
acoustic problems, such as the interaction of blade-tip shocks of a 
ducted, transonic rotor with wall boundary layers. 
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To solve equation (30/ subject to conditions (31) (with 
(31a) extended to the interval ■ • < X < 0), (32) and the Kutta 
condition at the trailing edge. In order that the solution be 
c ausal , i.e., that the displacement velocity vaves arise as a consequence 
of the notion of the shock (and are therefore present only for X ► - L) 
it is assuamd initially that o is appropriately situated in the upper 
half-plane to ensure that k j# k * both have positive imaginary parte. 

The final result is obtained by analytically continuing the solution on to 
the real w-axis. 

Set 


r* i(kX±y(k)x~) 

*1 " J C ± (k)e * • 


(A 


± according as *2 I 0 * Wiener-Hopf functional equations for € ± (k) 
are obtained by the usual procedure [22, Chapter 2], To simplify the 
argument condition (31a) can be replaced by 




X < 0 , % 2 - «■ 0 , 
(o • (3v ♦ KM) . 


(A 


Having obtained the solution ^(v), say, in this case, the cooplete 
solution is given by 


'i ■ i5 s * 0 r 


v 2 ( M -Uv)ftv<)) * ivt 


J__ V V,e iv 


(A 


Conditions (A2), (31b) lead to the functional equations 


Y(k)C * (k) 2 * (k-o+iO) 

“ L x (k) , 

(A< 

Y(k)C„(k) 

- L,(k) 

(A 



where X- 2 are regular in Im k < 0 end venieh es k ■* - i« • 
Conditions (32) respectively yield 
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c *< k) - c -< k > - W(w-lo) * V k > • <A6> 


T(k) fC/k) ♦ C.wj ♦ 2w(k~K -io) ■ °2 (k) • <"> 


where U^ U* ere regular in Zb k > 0 and vanish as k + ♦ i« . 
Using the nethod described by Noble [22] to eliminate 
one finds ultimately that 


Y(k)C 4 (k) 


1_ 1 
Ait k-o+xO 


* g /v 

V-r- T o 


+ 


Air 


.^(k-a+io) 


iA^^T 

k-K-iO 


(A8) 


The perturbation velocity and pressure will remain finite at the trailing 
edge (Kutta condition) provided y(k)C^(k) vanishes at least as fast as 
1/k as k ■*> • . This will be so if 


A - 


-x 


(A9) 


Inserting this into (A8) , substituting into (Al) and applying the integral 
operator (A3) gives the result (35) of the main text. 
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To solve equation (30) subject to condition (32) (vith 


B*0 


> i> 


3^2 

— - 0 for X < 0 „ % 2 - 0 , (A10) 

and the Rutta condition at the trailing edge* 

Set 

* 2 ■ *2 i * * 4 *2 • (All) 


where ** * 8 8i ven by equation (38). the *dge diffracted field 
when displacement velocity waves are ignored* First solve the reduced 
problem for , say, for which ** is replaced by 


20 


cos 


(X*2>' 


iYtt)X 


The total diffracted field is then given by 


(A12) 


*2 - f *2 Q WTO,)* iyWl dX . (A13) 

0 

Define 

n (* i{kX±y(k)x«) 

*2o " 4 I C < k > e * . (AU) 

i according as > 0 . This form automatically satisfies (32b) 

when B-0 . On * 2 ■ 0 , X < 0 equation (A10) implies that 


Y(k)C(k) - L(k) , 


(A15) 


where L(k) is regular in Imk < 0 and vanishes as k + - i- 


Condition (32a) givts 
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2C(k) + ini ' mtuL-m * 0<k) > <A16 > 

where U(k) la regular in Ink > 0 and vaniahea aa k + ♦ l*», • 

Solving (A15), (A16) by the Wiener-Hopf procedure and chooaing 
the value of A to aatiafy the Kutta condition, we find 


C(k) 


4wi^ 


|(k-K-JLO) 


- )] • 


(A17) 


Substitution into (A14) and application of the integral operator (A13) 
yields the result (40) of the audn text. 
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